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We present the key results from a comprehensive study of the refraction and focusing properties 
of a two-dimensional dodecagonal photonic "quasicrystal" (PQC), carried out via both full- wave 
numerical simulations and microwave measurements on a slab made of alumina rods inserted in a 
parallel-plate waveguide. We observe anomalous refraction and focusing in several frequency regions, 
confirming some recently published results. However, our interpretation, based on numerical and 
experimental evidence, differs substantially from the one in terms of "effective negative refractive- 
index" that was originally proposed. Instead, our study highlights the critical role played by short- 
range interactions associated with local order and symmetry. 

PACS numbers: 42.70.Qs, 41.20.Jb, 61.44.Br, 42.30.-d 
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Since the pioneering work by Yablonovitcbi and JohnS, 
"photonic crystals" (PCs) have elicited great attention 
from the scientific community, in view of the variety 
of peculiar electromagnetic (EM) bandgap, waveguid- 
ing/confinement, refraction, and emission effects attain- 
able through their use. Among the most intriguing appli- 
cations, it is worth mentioning those to negative refrac- 
tion and subwavelength imaging ( "superlensing" )^i^i^. 
The most typical PC configurations are based on dielec- 
tric inclusions (or voids) arranged according to periodic 
lattices in a host medium, and can thus be studied using 
well-established tools and concepts such as Bloch theo- 
rem, unit cell, Brillouin zone, equifrequency surfaces, etc. 

With specific reference to lensing applications, two dif- 
ferent approaches have been presented to obtain sub- 
wavelength resolution using a dielectric PC slab. In the 
first one, a PC with high dielectric contrast is tuned so 
as to behave (usually near a frequency band edge) like 
a homogeneous material with a negative refractive in- 
dex n = —1—, and the focus position of the flat lens 
follows a simple ray-optical construction^. In the sec- 
ond approach, "all angle negative refraction" (AANR) 
is achieved without an effective negative index, provided 
that the equifrequency surfaces (EFSs) of the PC are all 
convex and larger than the one pertaining to the host 
medium^. In this case, the focus position does not follow 
the ray-optical construction and is restrictec^. 

During the last decade, the discovery in solid- 
state physics of certain metallic alloys (the so-called 
"quasicrystals" ^"d^ ) whose X-ray diffraction spectra ex- 
hibit "noncrystallographic" rotational symmetries (e.g., 
5-fold or {K > 6)-fold, known to be incompatible with 
spatial periodicity) has generated a growing interest 
toward aperiodically- ordered geometries, leading to the 
study of the so-called "photonic quasicrystals" (PQCs). 
In this framework, useful tools for geometrical param- 
eterization can be borrowed from the theory of "aperi- 



odic tilings"™. Several recent numerical and experimen- 
tal studies have explored the EM properties of PQCs, 
in the form of two-dimensional (2-D) aperiodic arrays of 
cylindrical rods or holes, as well as 3-D structures fab- 
ricated via stereolithography (see [l^ and the references 
therein for a recent review of the subject). 

The study of PQCs entails significant complications, 
from both theoretical and computational viewpoints, as 
compared to standard (periodic) PCs. In spite of the 
lack of the aforementioned Bloch-type concepts and tools, 
approaches to the calculation of the density of states 
in PCQs have been proposed, relying, e.g., on rational 
approximantsi^ii^ or on extended zone schemes in the 
reciprocal spacers. However, many PQC properties and 
underlying mechanisms, which generally involve short- 
and long-range interactions as well as complex multiple 
scattering phenomena, are not yet fully understood. Nev- 
ertheless, results have revealed the possibility of obtain- 
ing similar properties as those exhibited by periodic PCs, 
with interesting potentials (e.g., richer bandgap struc- 
ture with lower and/or multiple frequencies of operation, 
higher isotropy, easier achievement of phase-matching 
conditions, etc.) in view of the additional degrees of free- 
dom typically available in aperiodically-ordered struc- 
tures. 

Recently, some examples of 2-D PQCs characterized by 
high-order (8-fold, 10-fold, and 12-fold) rotational sym- 
metries have been proposed as good candidates to exhibit 
negative refraction and subwavelength focusing effects 
with polarization-insensitive and non-near-field imaging 
capabilitie s^^d^ . Similar results were also obtained in 
case of acoustic wave&lS. Such effects and properties were 
interpreted in[l3 within the framework of "effective nega- 
tive refractive-index and evanescent wave amplification." 

In this letter, we present the key results from a compre- 
hensive study, based on full-wave numerical simulations 
backed by experimental verifications, which show that 
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FIG. 1: (color online) (a); Illustration of the Stampfli in- 
flation rule. Starting from the parent tiling represented by 
the gray-shaded central dodecagon, a big parent (red dashed 
lines) is generated by inflation, and filled up with copies of 
the original dodecagon placed at its vertices, (b): A portion 
of the tiling with the slab (~ 7a thickness) and wedge real- 
izations considered in[l3 (as well as in our study) marked by 
red and blue solid contours, respectively. The green-full-dot 
dodecagon corresponds to the parent tiling shown in (a), (c): 
Comparison (in a normalized scale) between the measured 
(black continuous curve) and simulated (red dashed curve) 
transmitted power for the PQC slab in (b) with a — 1.33cm. 



the "effective negative refractive-index" interpretation is 
questionable and unable to fully explain and predict the 
above effects, which instead arise from complex near-field 
scattering effects and short-range interactions critically 
associated to local symmetry points in the PQC. 

The PQC of interest is a conformally-scaled version of 
that in made of dielectric alumina rods (relative per- 
mittivity Cr = 8.6) of radius r = 0.4 cm placed (in air) at 
the vertices of a square-triangle 12-fold-symmetric tiling 
with lattice constant (tile sidelength) a = 1.33 cm, gen- 
erated according to the Stampfli recursive construction^^ 
(see also Fig. [T^) . Figure [TJa illustrates two particular 
slab- and wedge-shaped tihng cuts (as inll7|) considered 
in our study. 

In our numerical simulations, assuming the rods in- 
finitely long and parallel to the electric field, we employ 
a well-established full-wave technique (based on Bessel- 
Fourier multipolar expansion^^), which has been exten- 
sively applied to the study of 2-D finite-size PCs^^ and 
PQCs^. The experimental verification relies on mi- 
crowave (X-band) measurements on PCQ slabs made of 
alumina rods of height h = 1cm inserted in an aluminum 
parallel-plate waveguide terminated with microwave ab- 
sorbers, with a monopole antenna used to generate an 



FIG. 2: (color online) Simulated field intensity map at 8.836 
GHz for a PQC slab as in Fig. [TJ; (with lateral width of 42.9 
cm and thickness of 9.4 cm), illustrating the focusing of a 
source placed ai ya — and at a distance dx — 6cm from the 
slab surface. The dashed rectangle delimits the 10cm x 30cm 
area scanned in the image-side measurements. 



electric field parallel to the rods. The intensity /phase 
maps are collected using a HP8720C Vector Network An- 
alyzer and a computer-controlled x~y movable monopole 
antenna probe, in a setup similar to that described in 

As a preliminary check, we studied the transmission 
properties of the PQC slab in Fig. [Tb, for a fixed source 
position, within the 5 — 18GHz frequency range. Fig- 
ure [It compares the experimental and simulated (nor- 
malized) transmitted power through the slab. Three 
main bandgap regions are clearly observed, with fairly 
good agreement between simulations and measurements. 
We then studied the focusing properties within the fre- 
quency region 8—10 GHz (between the first and second 
bandgaps), by observing the field intensity maps at the 
image side. We found several frequency regions (includ- 
ing those reported infivl) where a clear focus was visible, 
with bandwidth ranging from few to thousands MHz. We 
note that this feature represents a first remarkable differ- 
ence with the (periodic) PC case, where the focusing re- 
gions tend to be more rare and well separated. The seem- 
ingly denser occurrence of focusing effects in PQCs could 
be attributed to their inherent self-similar natur^. We 
observed the most stable focusing effects around 8.350 
GHz. However, in what follows we concentrate on the 
focusing effects at 8.836 GHz, corresponding to the con- 
figuration considered in [TtI . In our study emphasis is not 
placed on quantitative assessments of the subwavelength 
focusing capabilities, but rather on the phenomenology 
interpretation, and its similarities and differences with 
respect to the periodic PC case. 

Figure [2] shows an example of a field intensity map in- 
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FIG. 3: (color online) As in Fig. (2] but details of the measured 
(a-e) and simulated (f-1) field intensity maps at the image 
side, for various source positions, (a), (f): Source at ys = 0, 
and at a distance dx = 6cm from the slab surface (cf. Fig. [2| ; 
(b), (g): = and = 7cm; (c), (h): ys ^ and d^ = 8cm; 
(d), (i): ys = 1cm and dx = 6cm; (e), (1): ys = 5cm and 
dx = 6cm. 



side and outside the PQC slab in Fig. [TJa. It is worth 
noticing that such slab is centered at the center of the 
tihng (x = y = 0, see Fig. [TJj), which is not only a center 
of local (12-fold) rotational symmetry, but also possesses 
reflection symmetries with respect to both x and y axes. 
The presence of a focus at the image side is clearly ob- 
served. However, from a comprehensive numerical para- 
metric study (see also 126' for more details), we found that 
neither a ray diagram approach can be used nor a prefer- 
ential propagation direction can be established to justify 
and predict the image formation. Conversely, we found 
that a local structure (in particular, the green-full-dot 
dodecagon parent tiling highlighted in Fig. [Ud) around 
the local symmetry center plays a key role. 

If the PQC slab behaved like a homogeneous mate- 
rial with a negative refractive index n = — 1, the source- 
image distance would remain constant and equal to twice 
the slab thickness^. Results for periodic PCs showing an 
almost isotropic refractive response (see, e.g., @) agree 
fairly well with this prediction. Moreover, in PC lenses, 
for lateral displacements of the source that preserve the 
distance from the lens interface as well as the structure 
periodicity, the focus remains unaffected, and its position 
follows the source location, in view of the absence of an 
optical axis. 

Figure [3] shows some representative measured and sim- 
ulated field intensity maps at the image side of the PQC 
slab, for orthogonal and parallel (to the slab interface) 
displacements of the source. Again, a general good agree- 
ment is observed between simulations and measurements. 
Specifically, Figs. [5^ and [3]: pertain to the configuration 




FIG. 4: (color online) Simulated intensity field maps at 8.836 
GHz for a coUimated Gaussian beam (with minimum spot- 
size ~ 2.8cm) normally impinging on several realizations of 
a PQC 60°-wedge. (a): Incident beam axis intersecting the 
local symmetry center, as in |i3[Fig. 2]. (b)-(d): Different 
wedge realizations obtained by random rigid translations of 
the tiling (so as to displace the local symmetry center from 
the incident beam axis). 



in Fig. [21 whereas Figs. and Figs. [3J:,[3Ji pertain to 

source displacements along the x— axis (i.e., orthogonal 
to the slab interface), which preserve the t/s = position 
(i.e., keeping the source facing the local symmetry center 
of the tiling). The focus position does not change sub- 
stantially. raising further concerns about the interpreta- 
tion in '13. Moving the source parallel to the slab surface, 
and therefore breaking the symmetry around the x— axis, 
one observes from Figs. [3jl,[3j and Figs. [3]3,[3j that even 
small displacements significantly affect the focus image, 
which undergoes a rapid deterioration until it completely 
disappears for the source placed at = 5cm. Interest- 
ingly, a focus can be still observed as long as one extends 
the slab size along the y— axis and places the source at 
ys — 13.67cm, i.e., directly facing the symmetry center of 
another big parent tiling (see FiglT]3). In this case, how- 
ever, the image (not shown for brevity) exhibits a worse 
quality. The complex interplay between local and global 
order and symmetry in the focusing properties, which 
can be glanced from the above results, is confirmed by 
parametric studies of PQC slabs of different thicknesses, 
where a variety of effects can be observed, ranging from 
localized (single and multiple) spots to beaming phenom- 
ena (see [2a for details). 

As a further check, we also carried out a numerical 
study involving a coUimated Gaussian beam impinging 
with normal incidence on the surface of PQC 60°-wedges 
extracted from the dodecagonal tiling. Results are shown 
in Fig. m Specifically, Fig. 2^ pertains to the wedge 
realization shown in Fig. [TJ), with the incident beam 
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FIG. 5: (color online) (a): As in Fig. [1 but for a PQC slab 
with lateral width reduced to 22cm. (b), (c): Details of the 
measured and simulated field intensity maps, respectively, at 
the image side. 



axis intersecting the local symmetry center. Again, our 
results are similar to those in [TtI . with the transmitted 
beam propagating mainly in the "negative" direction - 
a phenomenon previously interpreted within the frame- 
work of "effective negative refractive-index." However, a 
deeper study reveals that this effect too is critically re- 
lated to the mutual position of the incident beam and 
the local symmetry center. This is clearly visible in Figs. 
Sb'Stij pertaining to different PQC wedge samples ob- 
tained by random rigid translations of the tiling (so as 
to displace the local symmetry center from the incident 
beam axis), which display complex multi-beam features 
in the transmitted field, thereby highlighting the absence 
of a clear-cut refractive behavior. 

From the above results, which confirm the key role 
played by short range interactions involving a neighbor- 
hood of the parent tiling, one could speculate that the 
focusing properties of a PQC slab would be restricted 
to a limited range of incidence angles, and should there- 
fore occur also for significantly reduced lateral widths. 
Indeed, the focusing effects turn out to be rather robust 
with respect to lateral width reductions that do not affect 



significantly the modal field distribution in a neighboring 
region of the local symmetry center. Figure [5] shows the 
simulated and experimental results pertaining to a PQC 
slab with lateral width of only 22cm (i.e., six wave- 
lengths, nearly half of that in Fig. [2]), and yet still ex- 
hibiting a clear focus. Similar results were also obtained 
for PQC slabs with even smaller (only three wavelengths) 
lateral width, but with larger thickness^'. 

In conclusion, our numerical and experimental study 
of the refraction and focusing properties of dodecago- 
nal PQCs confirms some of the results reported in the 
recent literatureii, but shows that, contrary to the orig- 
inal interpretation, such results are not attributable to 
an "effective negative refractive-index." Instead, they 
arise from complex near-field scattering effects and short- 
range interactions critically associated to local symmetry 
points in the PQC, which were only glossed over in previ- 
ous studies. In this connection, it is worth recalling that 
local order and symmetry have already been observed 
to play a key role in a variety of PQC-related effects, 
including bandgap formation^'^, field localization^, and 
directive emission2§.. The new evidence and insights pre- 
sented in this letter constitute a further step toward a 
full understanding of the underlying mechanisms, which 
remains crucial for their judicious exploitation in the de- 
sign of novel compact optical systems. Within this frame- 
work, further numerical and experimental studies of var- 
ious PQC geometries are worth pursuing. 
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